Excitations of a Bose-Einstein condensate in a one-dimensional optical lattice 
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We investigate the low-lying excitations of a stack of weakly coupled two-dimensional Bose- 
Einstein condensates, which is formed by a one-dimensional optical lattice. In particular, we calcu- 
late the dispersion relations of the monopole and quadrupole modes, both for the ground state as 
well as for the case in which the system contains a vortex along the direction of the lasers creating 
the optical lattice. Our variational approach enables us to determine analytically the dispersion re- 
lations for an arbitrary number of atoms in every two-dimensional condensate and for an arbitrary 
momentum. We also discuss the feasibility of observing our results experimentally. 



I. INTRODUCTION 

O 

Cold bosons in an optical lattice provide a uniquely tunable environment to explore quantum phenomena. Some of 
, these phenomena have been known theoretically for quite some time, but with the advent of new experimental tools 
they have become a focus of attention. For example, Bloch oscillations of electrons in a metal are standard material 
in condensed-matter text books, but advances in the manipulation of cold atoms have made their experimental 
investigation also possible in this case Q. In addition, diffraction of matter waves by a pulsed optical lattice was 
studied by Ovchinnikov et al. Q. Optical lattices have also enabled the observation of some more exotic quantum 
phenomena such as number squeezing Q and collapses and revivals Q. Apart from these examples, Bose-Einstein 
condensates in optical lattices are particularly promising physical systems to study superfluid properties of Bose 
l— ~~ 1 , gases 0, |(| . Very importantly they realize the Bose-Hubbard model Q and can be used to investigate the quantum 
' phase transition from the superfluid into a Mott-insulator state 0, • This phase transition was recently indeed 
observed experimentally 0. 

■ In this paper we present a variational method to study the excitations in a stack of weakly-coupled two-dimensional 
(quasi)condensates. Such a system can be created by applying a relatively strong one-dimensional optical lattice to 

CO . an ordinary three-dimensional condensate. We focus on the transverse monopole and quadrupole-modes, but we also 
demonstrate how the method can be applied to study longitudinal excitations. We determine the eigenfrequencies 
of the monopole and quadrupole modes without any other approximations than those involved in our variational 
ansatz. In particular, this means that we can smoothly crossover from the noninteracting limit to the Thomas-Fermi 
regime. Moreover, the longitudinal wave length of the modulation is arbitrary, i.e., nearest neighbor sites can be 

■ completely out of phase. Using typical experimental parameters, we predict that the dispersion relations show a 
| strong dependence on the lattice potential. 

I i We consider the transverse excitations both for the ground state of the gas, as well as for the case that a vortex 
pierces through the center of each two-dimensional (quasi)condensate. We consider also the latter case, because a 
recent experiment observed the transverse vibrational modes of the vortex line in a trapped Bose-Einstein condensate. 
These modes are called Kelvin modes 10]. In the lattice three-dimensional effects, such as vortex line curvature, are 
expected to be less important. Therefore, we believe that Kelvin modes, and in particular their coupling to the 
transverse collective modes, are easier to study in a lattice than in a cigar-shaped three-dimensional condensate. Our 
results for a stationary vortex represent the first necessary step toward understanding the more complicated problem 
of the Kelvin modes of a vortex in an optical lattice. 

There is a number of important theoretical papers on the dynamics of Bose-Einstein condensates in optical lattices. 
For example, dynamical and modulational instabilities were studied in Refs. [Til fl^L fl3| and the adiabaticity of 
the nonlinear wave equations was explored by Band and Trippenbach |14|. Massignan and Modugno derived jl5| a 
relatively simple way to solve the three-dimensional Gross-Pitaevskii equation and used it to investigate the dynamics 
and expansion of the condensate in a one-dimensional optical lattice. Finally, Kramer et al. |16| used hydrodynamic 
equations to study the low-lying collective modes of a harmonically trapped Bose-Einstein condensate in the presence 
of a one- or two-dimensional optical lattice. In particular, they showed that the effect of the lattice is to renormalize the 
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interaction coupling constant and introduce an effective mass that accounts for the different inertia along the lattice 
potential. With these changes it is for instance possible to apply the results for harmonically trapped condensates 
obtained by Stringari |l7j . 

The study by Kramer et al. is somewhat related to ours. The most important difference, however, is that Kramer 
et al. are interested in different modes than we are. They deal with low-energy excitations along the long axis of 
the cigar-shaped condensate, i.e., the longitudinal modes. In the absence of a magnetic trap they correspond to 
the Bogoliubov modes with the familiar phonon spectra at large wave lengths. In the presence of a harmonic trap 
the spectrum becomes discrete and the lowest energetic modes are the center-of-mass mode and the (longitudinal) 
quadrupole-mode. The transverse modes we are interested in have superfluid flows orthogonal to the long axis of 
the condensate. Moreover, in the z-direction along the lasers of the optical lattice the condensate is for simplicity 
treated as completely periodic. As we see later, we take in first instance the atom number in each two-dimensional 
(quasi)condensate to be constant and equal in every site. This removes the Bogoliubov modes from the excitation 
spectrum, since they correspond to a density modulation propagating along the z-axis. If desired, they can however 
be easily incorporated into our approach. Indeed, at the end of the paper we briefly discuss the corresponding sound 
mode when the atom number in each two-dimensional (quasi)condensate is allowed to fluctuate. 

There are several experiments on condensates in a one- dimensional optical lattices Q, S S H E EU , but none 
of these address the problem we consider in this paper. Fort et al. ,2flj measured longitudinal excitation frequencies 
of the condensate in the presence of a one-dimensional optical lattice and this is the experimental paper most closely 
related to our work. In particular, for the breathing mode Fort et al. |20j do not report any dependence on the lattice 
depth and this is in agreement with our result as well as with the result of Kramer et al. • While we do not predict 
dependence on the lattice depth, we do predict dependence on the modulation of the excitation along the z-axis. In 
particular, we expect changes in the eigenfrequency if the sites are out of phase with each other and these changes 
can be large for typical experimental parameters. So far such phenomena have not been probed experimentally. 

The paper is organized as follows. In Sec. ^ we derive the theory used in this paper. In Sec. IIHI we apply this 
theory to a Bose-Einstein condensate without vortices and calculate the dispersion relation of the monopole and 
quadrupole- modes in the presence of a lattice. In Sec. IIVI we proceed by repeating similar calculation for the vortex 
state of the Bosc-Einstcin condensate. We end with a discussion of our results in Sec. [V] 



II. GROSS-PITAEVSKII THEORY 



Our starting point is a cigar-shaped Bose-Einstein condensate trapped by the potential 

7(r) = f KV+^V), (1) 

where u r and oj z are the radial and axial trapping frequencies, respectively, and M is the atomic mass. As we assume 
a cigar-shaped trap, we further have that u) z <C u> r . The condensate also experiences an one-dimensional optical 
lattice 

y (r)=U sin 2 (^), (2) 

where Vb is the lattice depth and A is the wave length of the laser-light. We assume that the lattice is deep enough 
so that it dominates over the magnetic trapping potential in the z-direction. When this is true and the number of 
lattice sites is large, i.e., A <C l z — \Z^/Muj z , we can in first instance ignore the magnetic trapping potential in the 
z-direction. 

The lattice potential splits the condensate into N s two-dimensional (quasi)condensates with a pancake shape. 
We assume that the lattice is sufficiently deep such that its depth is larger than the chemical potential of the two- 
dimensional (quasi) condensate |2l|. Using a Thomas- Fermi approximation for the two-dimensional (quasi)condensates 
we obtain a lower boundary that can be expressed as 



V » 2 9 / 7 



Na [jJ? 



4/7 



hu> r , (3) 



where N is the number of atoms per site and a is the three-dimensional scattering length. As a numerical example, 
we take an 87 Rb condensate in a trap with a radial trapping frequency Lu r /2ir — 100Hz and a laser-light wave length 
of A = 800nm. When the atom number in each site is between 100 and 1000, the lower bound on the trap depth Vo 
is between 0.05E r and 0.19-E r , where E r = ft 2 (2ir/\) 2 /2M is the recoil energy of an atom that absorbed one photon 
from the laser beam. 
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Although we are interested in a deep lattice, we consider here only the case that there is still full coherence across the 
condensate array. Specifically this means that the lattice potential should not be so deep as to induce a Mott-insulator 
transition. Typically the required lattice depth to reach the Mott-insulator transition in a three-dimensional lattice 
with a filling factor of one is of the order of 10E r . In a one-dimensional lattice the number of atoms in each lattice 
site is typically much larger than in a three-dimensional lattice and the transition into the insulating state requires a 
much deeper lattice. In mean-field theory the Mott-insulator transition in such a system occurs when Ur > 8NJ [2^ . 
where Ur and J are respectively the characteristic renormalized interaction and hopping parameters of the effective 
single-mode Bose-Hubbard model with Hamiltonian 



H = - J £ b]k 



<hj> 



2 ^ 



hi {h t - 1) 



(4) 



Using the same numerical values as in the previous paragraph, we estimate the critical lattice depth for the Mott- 
insulator transition to be between 56E r and 82E r , when the number of atoms in each site is again between 100 and 
1000. To the best of our knowledge the Mott-insulator transition in a one-dimensional optical lattice has not yet been 
observed. 

We use trap units from now on, i.e., the unit of energy is huj r , the unit of time is l/u> r , and the unit of length is 
l r = y/h/Mu> r . The Gross-Pitaevskii energy functional, which describes the system at low temperatures, is then 



^ 2 + v 2 ) 



Vo(r) , T 



2D 



TlUJr 



l*(r) 



|*(r)| s 



(5) 



where T 2B is the two-body T-matrix. In the above units the latter is related to the three-dimensional s-wave scattering 
length a through T 2B = A-Ka/l r . 

For a deep lattice potential it is natural to expand the condensate wave function in terms of wave functions that 
are well localized in the sites. More precisely, we expand 



* (r) = w ( z ~ z «) $ ™ ( x > V) ' 



(6) 



where n labels the lattice sites and z n — Xn/2l r is the position of the nth site. For now we do not specify the wave 
functions <&„ (x, y) of the two-dimensional (quasi)condensates, but for the wave function in the z-direction, w (z), we 
use the ground-state wave function of the harmonic approximation to the lattice potential near the lattice minimum. 
This harmonic trap has the frequency 



2vr 



ujl = —V2V /M 



and the wave function w(z) is thus given by 

w(z) 



exp 



211 



(7) 



(8) 



where II — y/h/MujL- 

Substituting the above ansatz into the energy functional and ignoring all but the nearest neighbor interactions, we 
get the energy functional 



d 2 r{ -^* n (x,y)\7 2 <S> n {x,y) + 



{ X 2 + y 2 ) + l ^\^n{x,yt 



\®n(x,y)\ 2 



J I d 2 r<$>* m (x,y)<$> n {x,y) } , 

<n,m> * 



where (n, m) indicates nearest neighbors, and 



U 2D =T 



2 LI 



dz\w{z)\ 1 ' 



7T / a 



2 \l 



(9) 



(10) 



is the two-dimensional coupling strength. Moreover, J is the strength of the Josephson coupling between neighboring 
sites and we have 



J = — J dzw*(z) 



2d^ 



w(z + \/2l r ). 



(11) 



4 



With these assumptions J is a time-independent experimentally defined parameter. Approximating the lattice poten- 
tial near its maximum by an upside-down parabolic potential we can calculate the Gaussian integral, with the result 



1 

8^2 



UJL_ 



IT 2 ' 



,-(V4ii) 2 



(12) 



The energy functional in Eq. @ is now almost two-dimensional. The third dimension is visible only in the last 
term that describes the coupling between neighboring layers. The energy is characterised by two parameters U2D and 
J, both of which are experimentally tunable. The importance of the on-site interaction term proportional to U2D 
can be enhanced by increasing the number of particles in the sites or by making the lattice deeper. Deepening the 
lattice also decreases the strength of the Josephson coupling J and makes the sites more independent. It should be 
noticed that while J is tunable, it is always positive. Physically this means that there is always an energetic penalty 
for having a phase difference between sites. 



III. EXCITATIONS OF THE CONDENSATE GROUND STATE 



In this section we study the transverse excitations of the ground-state of the stack of two-dimensional 
(quasi)condensates. Using a Gaussian ansatz for the wave functions of the (quasi)condensates, we solve the dis- 
persion relations for the monopole and the quadrupole-modes analytically. In Sec. IIII A1 we introduce the Gaussian 
ansatz and solve the excitations for an individial two-dimensional (quasi)condensate. In Sec. IIII Bl we proceed to 
calculate the band structure of the monopole and quadrupole-modes in the optical lattice. These sections also include 
technical details about the calculations. Such details are not repeated in Sec. IIVI when we consider the vortex state. 



A. Excitations for a single two-dimensional (quasi)condensate 

To account for the monopole and quadrupole modes of the two-dimensional (quasi)condensates in every site, we 
use a general Gaussian ansatz for the wave functions, i.e., 



$ n (x,y,t) = C„(i)exp 



(B xx<n (t)x 2 + B yy , n (t)y 2 + 2B xy ^ n (t)xy) 



(13) 



All three variational parameters Bij^ n (t) = B^ n (t) + iB'/j n {t) are complex. From now on we always use a prime to 
denote the real part of a complex quantity and a double prime to denote its imaginary part. The wave functions are 
normalized to the number of particles N at the site and therefore 



Cn{t) = \ ~ {BL,n(t)B' yy , n (t) - B' xy Jtf) 



2\l/4 



(14) 



As we fix the number of particles in every site, we are excluding the Bogoliubov modes propagating along the z-axis. 
It is, however, not difficult to account also for these modes as we show later on. The equations of motion for the 
variational parameters can be derived from the Lagrangian 



i = ! j d 2 r (5>;0r, y ,i) 

\ n 



d$ n (x,y,t) d$* n (x,y,t) 



at 



dt 



(15) 



Let us first investigate the behaviour of an individual two-dimensional (quasi)condensate. Without the interlayer 
coupling the part of the Lagrangian quadratic in the deviations (t) turns out to be equal to 



1 



U 



_ = _ ( f > i" + f ' f" + Of' e" ) + — 

at>2 V-xx^xx ~ ^yy^yy ' ^xy^xyj ~ jj 



45, 



Bo 



' xy 
2 



("253 
1 



— + — ^ (V 2 + e' 2 + e' e' + e' 2 ) + ( — + —] (e' + e' Y 
g3 ~ 2BqJ yy yy y \4Bg 4-Bq/ 



(e' +e' "| 2 + 2e' e' + e " 2 + e" 2 + 2e" 2 

\^xx 1 yy } 1 xx yy 1 ^xx 1 yy 1 xy 



(16) 
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where 

V27T \lrj V OJr 

We also defined the equilibrium solution of Bn(t) as B$. Hence, Bij(t) — Bo&ij + (t). We also suppressed the site 
index n. The equilibrium solution Bq is given by minimizing the zeroth-order term of the energy, i.e., 

Eo = ^ + ^ + UBq (18) 

with the result 



In Eq. (|16|l we show only the part relevant for the dynamics and we ignored the zeroth-order term whose minimisation 
lead to the result in Eq. fTt^l . 

We are now in a position to find the frequencies for the collective excitations we are interested in. Let us start with 
the monopole to = mode, which is alternatively also called the breathing mode. For the monopole-mode we can set 
£xx — tyy = e and e xy — 0. With this choice the Lagrangian is greatly simplified to 



L _ 1 
N ~ 2B~l 



^ J' 2 R <:" 2 



Bo 



(20) 



The equations of motion for e' and e" are the Euler-Lagrange equations that result in two first-order differential 
equations 

(' + 2B e" = 

-e» + A e ' = . (21) 
Bo 

These equations can be cast into a single second-order differential equation for e' 

e' = -4e', (22) 

which describes sinusoidal oscillation with a frequency 2. The frequency of the monopole-mode is therefore 

loo = 2 (23) 

and it is independent of the strength of interactions. This is in agreement with previous results |23j . 

The quadrupole to = ±2 modes are captured by the choice e xx — —tyy — £• We then have just two (complex) 
variational parameters, e and e xy . In the Lagrangian in Eq. (|16fl there are no terms that couple e to e xy . Therefore, 
the dynamics of these parameters separates and both turn out to have the same oscillation frequency. Above we gave 
the necessary technical details in the derivation of the monopole-mode frequency. As the quadrupole-mode frequency 
can be dealt with in a similar fashion, we simply give the result. The quadrupole-mode frequencies are given by 



lo±2 = ^/2 + 2BI (24) 

For the ideal-gas Bo = 1 and the quadrupole frequency is again 2. In the Thomas- Fermi limit Bo tends to zero and 
the quadrupole frequencies app roach lu±2 — ■> y/2. Again this result is as expected (171 124 l25T | . Our treatment also 
captures the scissors mode [2Gl . but in the axial symmetric case we are considering here the scissors mode turns out 
to be degenerate with the quadrupole-mode. 

Incidentally, it should be remembered that the degeneracy of the quadrupole modes is lifted in a rotating trap. If 
the trap is rotating with frequency f2 around the z axis we should include a term —Q,{L Z ) into the energy functional, 
where (L z ) is the expectation value of the angular momentum component in the z-direction. The angular momentum 
of the equilibrium solution is zero and the new term will only contribute in second order. The new contribution to 
the energy is 

-n{L z ) = ^[e' xy e"-el y e']. (25) 
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This term couples the dynamics of e and e xy , but the resulting 2x2 matrix problem is easy to solve. The quadrupole- 
mode frequencies in a rotating trap are 



oj± 2 = V2 (l + Bl) 1/2 ±V2Q 



(26) 



From this result it is clear that the quadrupole mode with m — — 2 becomes thermodynamically unstable when 
il > + Bq) /2. This result corresponds to the Landau-criterion for the quadrupole-modes, and has been shown 
to play an important role in the nucleation of vortices [53, H^, HH H3, []n| ■ 

B. Influence of the lattice on the excitation frequencies 

We are now in the position to discuss the influence of the lattice potential. To make progress we must determine 
the coupling integral 

Imn = J d 2 r<S>* m { Xl y)<S> n {x,y) (27) 
to a sufficient accuracy. This will contribute to the energy a Josephson coupling 

Hj = - J ^ Knn: ( 28 ) 
(n,m) 

where the (n, m) indicates nearest neighbors. Here the imaginary part of I mn is not relevant since its contribution 
to the energy vanishes when the sum over the nearest-neighbors is calculated. For the monopole-mode we get up to 
second order in the deviations the result 

V 1 11 

~ 1 g^>2 Vn 1" fc m "T" At n T zt m / ' ^2 t n t m 1" ^J^- n m ^ ' 

and for the quadrupole-mode we have 

^ = 1_ 8^ (l e ™| 2 + l e ™| 2 + l e ^™| 2 + l e ^™| 2 ) + 

4- 1 (V F ' 4- f"f" 4- f' f' 4- f" f" 1 fSM 

1 V n m 1 L n u m 1 xy,n xy,m 1 ^xy,n^xy,mj ' 

In these formulae the first subindex of e^.n identifies the variational parameter in question and the second indicates 
the lattice site. For identical nearest-neighbor wave functions the overlap integral I mn should be exactly TV, which is 
indeed the case in both equations (|29|) and 1|3U|) . 

Some terms in Eqs. (|2*§1) and i|3U|) are purely on-site, but terms of the type e„e m are not. This complication is 
remedied by going to Fourier-space. We define the Fourier transform in such a way that the function /„ in coordinate 
space is expressed in terms of its transform fk as 

fn = -7== ^2 exp[ikz n ]f k . (31) 

Here iV s is the number of lattice sites which we, for notational convenience, assume to be an odd number. Moreover, 
k is the wave number and the lattice spacing is d — A/2. 

First we transform the diagonal terms in the Lagrangian. For example 

E % = E w E m« ex p ^« ( fc + fc ')i = E = E i^i 2 < ( 32 ) 

n n s k,k' k k 

where the sum over the lattices sites n gave the Kronecker delta 5k\—k which removed one of the momentum sums. 
The last step is a result of the fact that f n was a real function, so ft = f-k- Nearest-neighbor terms are somewhat 
more complicated. As an example 

fnfm = gTf E E -f k f k ' t cxp ( fcz ™ + k ' z n+i)) + cxp (i (kz n + k'z n -i))] . (33) 

<n,m> s n k,k' 
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We can perform the sum over n and get 

]T Mm = cos ( fc ' A / 2 ) fkfk>h>,- k = J2 cos ( fcA / 2 ) I A 



(34) 



In Fourier-space the Josephson coupling generally thus introduces factors of cos(fcA/2) — 1 into the Lagrangian. 

Now that we know how to transform to Fourier-space, we can proceed to derive equations of motion for each value 
of the wave vector fc. Since two different values of the wave vectors do not couple, this is not technically any more 
complicated than our previous treatment of an individual (quasi)condensate. The equations for each wave vector can 
be solved separately. We demonstrate this again for the simplest case, namely the breathing mode. Let the Fourier 
transform of e n be In Fourier space the Lagrangian for the breathing mode is 



1 - -^2iE e fe 



B — 2 J [ cos 



kX 



"|2 



\4\ 



(35) 



Keeping in mind that e_fc = t* k we get equations of motion for e' k and e'L For example by considering the variation of 
the Lagrangian with respect to e" k we get 



4-2 



Bq — 2 J cos 



kX 



e k = 



(36) 



and by considering variations with respect to e'_ k we get the differential equation for e' k 



-k ~ J ( cos (t 



0. 



(37) 



The dispersion relation for the monopole-mode can now be simply read out from this pair of equations. The 
quadrupole- modes can be dealt with in the same way although the equations are somewhat longer. 

For convenience we assume that the contribution from terms proportional to J 2 are very small. With this simplifi- 
cation we get the dispersion relations for the monopole and quadrupole-modes 



w (fc) = 2 



i / ( Bo + — 



kX 

cos — 1 

2 



1/2 



(38) 



LU ±2 {k) = V2 



1 + BZ - J [ 3B + — 

-DO 



cos 



"I 1/2 



(39) 



We emphasize that our results where terms under the square root proportional to J 2 are ignored should be used 
with some caution. The terms proportional to J 2 are not always negligible compared to the other contributions. In 
particular if the trap depth or the on-site number of particles is small, there is a range of experimentally relevant 
parameter values where terms proportional to J 2 can be relatively large and should be included. They will not change 
the qualitative behavior of the dispersion relations, but can affect quantitative results. While we choose to work in 
the regime where terms proportional to J 2 are small, it is not difficult to include these missing terms. For example, 
Eqs. (|36|l and (|37|l show that the exact frequency for the monopole-mode obeys 



wg(*) = 4 



Br. - 2 J I COS 1 — 1-1 



(40) 



In Fig. ^ we show the dispersion relation for the monopole-mode as a function of k and U. 

In the limit of long wavelengths it is permissible to expand the cosine factors. The excitation has then the same 
dispersion as that of a free-particle A + h 2 k 2 /2m* with some effective mass m* and a gap A. For the monopole-mode 
we, therefore, predict an effective mass 



4B 
J(B 2 +2) 



,A 2 



(41) 
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FIG. 1: monopole-mode frequency as a function of k and U, when J — 0.05. The surface in this figure was calculated using 
Eq. l$g) . 



and for the quadrupole-mode we get 



4^0 y/2 . , , 
m±2= J(3Bg + l) l^J' (42) 



It is quite interesting to observe that the effective masses of different modes are different. In particular, the effective 
mass of the quadrupole-mode is bigger than the effective mass of the monopole-mode. This can be understood by 
considering the overlap integral between neighboring sites. For the monopole-mode the coupling between the nearest 
neighbors is to a large extent determined by the integral 

d 2 r\<t> (x,y)\ 2 (x 2 +y 2 ) 2 , 

where <&o( x il/) is the equilibrium wave function. In contrast, for the quadrupole-mode the coupling is determined by 
the integral 

' d 2 r\<P (x,y)\ 2 (x 2 -y 2 ) 2 . 

It is clear that the latter integral is smaller than the first one. As the effective mass is inversily proportional to the 
strength of the nearest-neightbour coupling, the quadrupole-mode therefore has a bigger effective mass. 

The fact that the dispersion relation is similar to the free-particle dispersion relation is reflected in the dynamics. 
A sinusoidal modulation moves with velocity v ~ h(k)/m* in the z-direction. In a finite system a pure sinusoidal 
modulation is not possible and the excitation corresponds to a wavepacket centered around (k) and with some nonzero 
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width Ak. If the system is large enough, i.e., much bigger than 2n/Ak, the width of the packet can be small and one 
should be able to observe such propagation before the excitation hits the outer edge of the condensate. 

More generally we can expand the dispersion relation around any value of the wave vector. In terms of a function 
C{J,Bq) that depends on the mode in question, the excitation energy up to lowest-order in J looks like 



u(k) = lu(0) + C (J, B ) 



Expanding this expression around ko we get 
uj{k) w w(k ) — C (J, B Q ) 



A . / k X\ A 2 (k X \ 2 

sin | — — I [k — ko) + — cos — (k - fc ) 



(43) 



(44) 



When k = we get the results for the effective masses we presented earlier, but some special cases are also of interest. 
In particular, when kX/2 = n we obtain the same expansion as with k = 0, but the constant in front of (k — ko) 2 has 
a negative sign. In this regime the effective mass is therefore negative. In the regime of a negative effective mass one 
encounters modulational instabilities as discussed in Refs. [TTL IT3. fT3| . 



IV. EXCITATIONS OF THE VORTEX STATE 



In this section we consider a system of weakly-coupled two-dimensional (quasi)condensates that has a vortex piercing 
through the center of each (quasi)condensate. For such a system our earlier ansatz in Eq. (|13|) is inadequate. For an 
individual condensate it is known that the presence of a vortex should not change the dispersion of the monopole- 
mode, but it will lift the degeneracy of the quadrupole-modes. Physically this is due to the fact that the quadrupole 
excitation, depending on the sign of the quantum number m, travels either in the same direction of the superfluid 
flow or opposite to it. As the monopole-mode is easier to tackle than the quadrupole-modes, we start with that in 
Sec. IIV A1 In Sec. IIV 51 we solve for the quadrupole-modes of an individual two-dimensional (quasi)condensate and 
in Sec. IIV O we include also for the quadrupole-mode the optical lattice into our discussion. 



A. monopole-mode in the presence of a vortex 



The vortex state has a superfluid flow around the vortex core. This flow diverges in the core and for this reason the 
density of the condensate must vanish in the vortex core. The simplest ansatz having these two desired properties is 
(in polar coordinates) 



$ n (r, 4>) oc rexp [i<fi] exp 







e n (t)r 2 ' 


2 


exp 


2 



(45) 



The ansatz is almost the same as in the previous section for the monopole-mode of the state without a vortex. The 
only differences are the first two factors that give the vortex the properties we were after. The size of the vortex core 
region in Eq. (|45(l is about 1/^/Bq and it does not diminish as the number of particles is increased. This is in principle 
incorrect, since the length scale for the vortex core size is set by the coherence length and the coherence length in the 
center of the condensate gets smaller as the number of particles is increased. We expect that this unphysical behaviour 
close to the vortex core is not relevant to the physics of the collective modes at hand. In the end of the calculations we 
can reproduce the known results for the individual pancake to a good accuracy and thus our expectations are indeed 
well justified. 

Using similar techniques as for the condensate without a vortex, we calculate the monopole-mode of a condensate 
with a vortex as 



Lu (k) 



I :-!•/ ! Bo + ^~ 

t>0 



k\ , 

T 1 



n 1/2 



(46) 



where we have again assumed that J 2 terms under the square root can be ignored. The equilibrium solution is now 
given by 



Bo 



U 



(47) 
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The result is similar to Eq. (|38(l . but the constant in front of the cosine term is different indicating a difference in the 
effective mass. For the monopole-mode in the presence of a vortex we get 



///,, 



45 



iJ{Bl + l) \uj r X 



(48) 



We can see that the effective mass of the monopole-mode of the vortex state is somewhat smaller than the effective 
mass in the absence of a vortex. This can be understood by comparing the relevant overlap integrals for the wave 
functions with and without the vortex. Since the (quasi)condensate wave function with a vortex is more extended 
than without a vortex, the strength of the nearest-neighbor coupling is increased and, therefore, the effective mass is 
reduced. 

B. quadrupole- modes of the single pancake in the presence of a vortex 

As we mentioned before the quadrupole-modes are more complicated. For the quadrupole-modes we use the ansatz 



$„ (r, <fi) oc r exp [i<f>] exp 

~ r exp [i(j>] exp 
- sin (20) ( 



B r 2 



2 

B r 2 



exp 



1 + cos 



[1 + a exp [—2iq 



m ( 



a r 

2 



, Cx V 2 



(49) 



where a denotes a new variational parameter and the last expression is an expansion of the first line to the first order 
in the deviations. This ansatz looks somewhat complicated, but this is needed to build in the relevant physics. This 
is most easily seen by considering the noninteracting limit where the wave functions are known analytically. 

In the noninteracting limit the vortex states with angular momentum projections equal to ±N are degenerate. This 
implies that linear superpositions of these states have the same energy. As a result, there exists a quadrupole-mode 
with zero frequency in this limit. To capture this mode the variational parameter a is included in the ansatz. To 
understand this, assume that a — and expand the exponent in Eq. We get 



exp 



(x 2 - y 2 ) 



e r 



cos (2q 



■ sin (20) + O 



(e 2 e 2 



CC 



xy) 



(50) 



For clarity assume also that e is real and e xy = —ie. The disturbance then couples to the wave function 



r exp [i(j>] exp 



B r 2 



r 2 exp [— 2i<i 



which is the wave function of the anti- vortex state multiplied by r 2 . This state has obviously a different energy than 
the true anti-vortex wave function. As a result the ansatz without a gives a wrong frequency for this mode in the 
ideal-gas limit. To avoid this problem we need the additional variational parameter a to give a nonvanishing amplitude 
for the correct anti- vortex wave function in the ideal-gas limit. The fact that this ansatz really couples to the correct 
anti- vortex wave function is most clearly seen by setting e = e xy = in Eq. (|49|l . Otherwise the ansatz is very similar 
to the ansatz we used for the (quasi)condensate without the vortex. In the noninteracting limit the m = +2 mode 
requires coupling to a wave function with angular momentum m = 3 and with a small distance behaviour that should 
be proportional to r 3 . In Eq. (|49|l this is indeed the case, as can be verified by setting a equal to 0. 

Using the above ansatz we can determine the quadrupole-modes of a single pancake analytically for the full pa- 
rameter regime from the noninteracting limit to the Thomas-Fermi regime. The equilibrium solution is the same as 
for the monopole-mode and is given by Eq. (|47|l . To second order in the deviations the various contributions to the 
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Lagrangian are 



L T = -2a'a" - i (e'e" + 4,4.) + ^ [«' (*" + i'„) + *" (?„ - f 



1 / i i O i i O . ii . ii i \ 1 



Lv = 7753 (l £ | 2 + M') - o~r~ K ( e ' - e "v) + a " ( e " + e -y)] 



Lk = 4^ (|e|2 + |e * y|2 + 4e ^ " 4e " 4 ^ + 2 [a ' (e ' _ + a " (e " + (51) 

i« = C/Bo {hi 2 + g|j ~ e" 2 + 4 - 4' 2 ) - ^ [«' (5e' + e%) + a" (5e' xy e")) J , 

where is due to the part of the Lagrangian containing the time-derivatives, Ly is due to the potential energy, Lk 
is due to the kinetic energy, and Ljvl is the contribution due to the interactions between atoms. 

With this result we can solve for the eigenmodes of the system. The problem is essentially that of solving the 
eigenvalues of a 3 x 3 matrix. This matrix has three (generally) nondegenerate eigenvalues and two of these correspond 
to the quadrupole-modes. The third mode is of no interest to us here. In this modes the deviation from the equilibrium 
is a superposition of various trap states, among which the m = 5 component has an incorrect short distance behaviour 
that causes the energy of this mode to strongly increase with increasing atom number. The frequencies of the 
quadrupole-modes can be calculated analytically, but the results are too long to be given here. However, they do 
not cause any computational problems. In Fig. [21 we show the frequencies of the quadrupole-modes based on our 
ansatz and compare them against the values computed numerically with the Bogoliubov-deGennes equations |32| . 
The agreement is very good over the whole range of interaction strengths. 

In the limit of a nearly ideal-gas the quadrupole frequencies are w_2 = U and ll>2 = 2. For large atom numbers the 
quadrupole-mode frequencies are given by 



^±2 



= V2±-i= (52) 



and the splitting between the modes is U2 — — y/2/U. Zambelli and Stringari [33j | used sum rules to show that 
the splitting between the quadrupole-modes in the limit of large atom numbers should be 

UJ 2 - UJ-2 = / 2 \ • ( 53 ) 

Here L z is the z-component of the angular momentum operator. With our ansatz we have (L z ) = 1 and (r 2 ) = 2/Bq 
and the splitting of the quadrupole-modes is indeed the same as the result based on the sum rule approach. 



C. Influence of the lattice on the quadrupole-modes of the vortex state 

Including the lattice structure makes the already complicated equations even more complicated |34||. In Fourier 
space the nearest-neighbor interaction introduces a new term to the Hamiltonian 




where a' k is the Fourier transform of the real part of a n , e' k is the Fourier-transform of the real part of e, and e' k 
is the Fourier-transform of the real part of e xy . Similar notation applies to the imaginary parts and for example 
|efe| 2 = e' k e'_ k + e' k ('L k - In Figs [j^and Q]we show in detail the resulting dispersion relations for the quadrupole-modes 
as a function of k and interaction strength U . 

Even though the general formulae are too complictated to be given here, the ideal-gas limit and the Thomas- 
Fermi limit give us simple formulae. In the limit of weak interactions we have cj_2 = U — 2J (cos (fcA/2) — 1) and 
u>2 = 2 — 2 J (cos (fcA/2) — 1) and in the limit of strong interactions or large particle numbers we have 

W±2,d = V2 
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FIG. 2: Splitting of the quadrupole-modes for the Bose-Einstein condensate with a vortex as a function of interaction strength. 
The solid line is the analytical result based on the ansatz in Eq. 1)49^ and the open circles are calculated by solving the 
Bogoliubov-deGennes equations numerically. 



In the limit of large particle numbers the effective mass of the quadrupole-modes thus becomes 

16 / h 



l ±2,v 



3JVU \u r \ 



(56) 



This result indicates that we expect the quadrupole-modes to have about three times larger effective mass than the 
breathing mode. Again this can be understood by overlap arguments. 



SUMMARY AND CONCLUSIONS 



We have calculated the band structure of the most important transverse collective excitations of a stack of two- 
dimensional Bose-Einstein condensates in a one-dimensional optical lattice with and without a vortex. Our variational 
approach enables us to crossover smoothly from the ideal-gas into the Thomas-Fermi regime and to treat the interlayer 
coupling without other approximations that those involved in the variational ansatz. We have also calculated the short 
wave length part of the excitation spectra, which means that in our approach neighboring sites can be completely out 
of phase with each other. Using our general results for the excitation frequencies, we derived predictions for the the 
effective mass of the monopole and quadrupole-modes. We noticed that the effective mass is sensitive to the mode 
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FIG. 3: Dispersion of the quadrupole-mode with m — —2 for the Bose-Einstein condensate with a vortex when J = 0.05. The 
figure is based on the wave function ansatz in Eq. 1491 . 



in question as well as to the presence of a vortex. In this paper we have only focused on the linear response of the 
system. For large modulations nonlinear effects can become important |llill2l ll3j. In particular, assumptions about 
a nearly homogeneous condensate can break down as the system becomes dynamically unstable towards large density 
modulations. 

Experimentally the kind of excitations we have discussed in this paper can be created by modulating the radial 
trapping frequency u r as a function of z. One possible way to excite the monopole modes is to have two counter- 
propagating laser beams with a Gaussian intensity profile. Due to the optical dipole force the intensity profile of each 
one of the beams would provide the trapping in the radial direction, while the interference between the beams would 
provide the necessary modulation. To excite equal superposition of m = ±2 quadrupole-modes sheets of laser-light 
can be considered. 

In the limit of large interactions the constant in front of the k dependent part of the dispersion relations always 
scales with J\f~U . This number is a good measure of how strong the effects due to the lattice are. If this number is 
small, the lattice effects arc hard to distinguish experimentally from the dominant single site result. In Fig. [5]we plot 
JVU as a function of the on-site atom number for a few different lattice depths. As can be seen, the effects of the 
lattice for the modes we are considering can be very pronounced and should be easily observable. 

In the system we have discussed in this paper, each (quasi)condensate becomes very quickly two-dimensional 
as the depth of the lattice is increased. In particular the coherence length in the center of the two-dimensional 
(quasi)condensate quickly becomes larger than the thickness of the pancake. In low dimensions phase fluctuations 
are expected to be more pronounced |35i ISfiL l3^| . In our treatment we ignore such fluctuations. In a two-dimensional 
system there is a true condensate at zero temperature and then the phase fluctuations are not expected to play 
a major role. At nonzero temperatures phase fluctuations become more important, but are expected to be more 
pronounced between sites that are well separated. In our parameter regime the tunneling term coupling neighboring 
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FIG. 4: Dispersion of the quadrupole-mode with m — 2 for the Bose-Einstein condensate with a vortex when J = 0.05. The 
figure is based on the wave function ansatz in Eq. 1149^ . Note that for clarity the viewing angle is different from the previous 
figures. 

sites will establish phase coherence between neighbors. As a result the two-dimensional (quasi)condensates are not 
strictly two-dimensional since they "see" the third direction through the tunneling term. As the distance between the 
sites increases the phases become less correlated, but as we are only interested in the nearest-neighbor couplings, such 
effects are not important. Consequently, we expect our model to be applicable also at small but nonzero temperatures. 
Phase fluctuations may cause a slight reduction in the strength of the Josephson coupling, but would leave our results 
otherwise unchanged. 

In this paper we have choosen to fix the number of atoms in every two-dimensional (quasi)condensatc. In our 
variational approach it is not difficult to include atom number as well as global phase fluctuations by replacing in 
our variational ansatz y/~N by -J 'N n {t)e Wn W , where N n denotes the number of atoms and v n the global phase of the 
(quasi)condensate in every site. In the simplest case where we neglect the couplings with the transverse modes, we 
find that at long wave lengths there exists a phonon mode with the sound velocity 




which agrees exactly with the results obtained previously d, 0, ll^l • 

In a recent experiment the Kelvin modes of a Bose-Einstein condensate with a vortex were observed [Tjj- In the 
model that we have presented in this paper the vortex is always in the center of each pancake. In the future we plan 
to relax this condition and consider also the dynamics of the vortex. In this manner it is possible to study the Kelvin 
modes in an optical lattice, and in particular their coupling to the transverse excitations which were our main focus 
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FIG. 5: The quantity JVU as a function of the number of 87 Rb atoms in the each lattice site for different lattice depth. The 
solid line is for the depth Vb = SE r , the dashed line is for Vb = 9E r , and the dot dashed line is for Vb = WE r . The wave length 
of the laser-light was taken to be A = 800nm and the radial trapping frequency was L0 r /2-K = 100Hz. 



here. 
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